Self-induced and induced transparencies of two-dimensional and 

three- dimensional superlattices 



Yuriy A. Romanov and Julia Yu. Romanova 
Institute for Physics of Microstructures RAS, 603600 Nizhny Novgorod, 

Russia 

Lev G. Mourokh and Norman J.M. Horing 
Department of Physics and Engineering Physics, 
Stevens Institute of Technology, Hoboken, NJ 07030 
(February 1, 2008) 

Abstract 

The phenomenon of transparency in two-dimensional and three-dimensional 
superlattices is analyzed on the basis of the Boltzmann equation with a colli- 
sion term encompassing three distinct scattering mechanisms (elastic, inelas- 
tic and electron-electron) in terms of three corresponding distinct relaxation 
times. On this basis, we show that electron heating in the plane perpendicular 
to the current direction drastically changes the conditions for the occurrence 
of self-induced transparency in the superlattice. In particular, it leads to 
an additional modulation of the current amplitudes excited by an applied 
biharmonic electric field with harmonic components polarized in orthogonal 
directions. Furthermore, we show that self-induced transparency and dynamic 
localization are different phenomena with different physical origins, displaced 
in time from each other, and, in general, they arise at different electric fields. 
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I. Introduction 



Semiconductor superlattices have been at the focus of attention for several decades, due 
to their unique electronic properties. The additional spatial periodicity of the superlattice 
leads to the formation of narrow Brillouin minizones and energy minibands [JJ-|5]. Bloch 
oscillations || and Wannier-Stark levels can be observed in superlattices due to the nar- 
rowness of these minibands even in relatively weak static electric fields (10 2 — 10 A V/cm). 
The Bloch oscillations are due to Bragg reflections by the periodic superlattice potential and 
are characterized by the frequency Q c = eE c d/h and amplitude Z c = A/2eE c , where E c is 
the constant electric field applied along the axis of the superlattice of period d and miniband 
width A. In the case of an applied harmonic ac field, the Bragg reflections do not generate 
a new type of oscillation beyond that of the static field, but they do modulate electron 
motion during the field period. This modulation is described by oscillatory dependencies of 
the amplitudes of electron velocity harmonics on the amplitude, E\, and/or the frequency, 
Ui, of the applied harmonic field 0,0] . Manifestations of this modulation can be found in 
various nonlinear macroscopic effects and, in particular, in superlattice transparency []7|-|TT||. 
The zero th harmonic of electron nonlinear oscillations responsible for dc current is of special 
interest. Its vanishing corresponds to electron spatial localization and is called dynamic 
localization |J. This dynamic localization occurs only for electrons having a sine-like dis- 
persion law and for specific ratios of amplitude and frequency of the applied field, such that 
J (eEid/TiUi) = (where Jq(x) is the zero th order Bessel function). In the case of devia- 
tion from a sinusoidal dispersion law, dynamic localization can arise only at multifrequency 
fields [0. In the literature (see, for example, Ref. fl2]|), dynamic localization is very often 
erroneously identified with self-induced superlattice transparency, predicted in Ref. || and 



verified experimentally in Ref. [10]. It was shown there that the macroscopic polarization 
of the electron gas can vanish and the superlattice behaves almost like a linear dielectric 
having the permittivity of the crystal lattice in the absence of mobile electrons, with small 
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nonlinear absorption due to electron oscillations. The conditions for this effect within the 
single r-approximation for a one-dimensional superlattice sample are the same as for dy- 
namic localization, but they have different physical origins. The error of their identification 
was shown in Refs. |7|JT3[] for a three-dimensional sample with a one-dimensional superlattice 
and it will be further confirmed below for two- and three-dimensional superlattices. 

The single r-approximation has been reasonably successful in describing the cases of a 
one-dimensional superlattice and also a one-dimensional model of a three-dimensional super- 
lattice without redistribution of energy and momentum among the various degrees of freedom 
due to electron scattering. However, as was shown in Ref. fljjjlsf , such energy-momentum 
redistribution resulting from scattering can substantially affect superlattice properties and, 
in particular, the current-voltage characteristics can change due to transverse heating. To 
overcome the deficiency inherent in the single r-approximation, we develop a new method 
based on the Boltzmann equation with a collision term encompassing three distinct relax- 
ation times. The three relaxation times include (a) a time for redistribution of energy and 
momentum supplied by an electric field to a given electron among its various degrees of free- 
dom, (b) a time for redistribution of energy and momentum among all electrons by inelastic 
electron-electron interactions, and (c) a time for transfer of the excess energy to the crystal 
lattice. In this, we employ a separation of the relaxation processes into elastic, inelastic 
and electron-electron, which is commonly recognized in the study of nonlinear properties of 
semiconductors at high fields (see, for example, Ref. fl6||PT| ). The resulting balance equa- 



tions which we obtain can be solved analytically for systems having high symmetry (two- 
and three-dimensional superlattices). However, the qualitative results obtained here are 
also valid for bulk semiconductors having a one- dimensional superlattice, i.e. for structures 
with minibands in the growth direction and free electron motion in the lateral plane, which 
can be studied experimentally currently. Moreover, at the present time, three-dimensional 



cluster lattices are actually grown []n| and technological progress p9|-f2T[] offers hope that 
two-dimensional and three-dimensional superlattices will be fabricated in the near future 
using quantum dots, relating directly to our studies here. Furthermore, a simple three- 



dimensional Kronig-Penney model was proposed in order to describe such quantum dots 
superlattices theoretically |[22|| . On the basis of the three relaxation time description, we 
show that electron heating in the plane perpendicular to the current drastically changes 
the self-induced transparency of the superlattice. In particular, it leads to an additional 
modulation of the current amplitudes excited by an applied biharmonic electric field with 
harmonic components polarized in orthogonal directions. We obtain analytical results in the 
weak scattering approximation (lot » 1) and extend the analysis numerically for stronger 
scattering. 

This article is structured as follows. In Section II, starting from a Boltzmann equation, 
we derive balance equations for average electron velocity (current) and electron energy by 
means of the new collision term accounting for momentum and energy redistribution among 
the various degrees of freedom. On the basis of these balance equations the theory of self- 
induced transparency and its mechanisms is developed in Section III for one-, two- and 
three-dimensional superlattices in the presence of a high-frequency (hf) harmonic field. In 
Section IV we analyze the amplitude modulation of hf current by an orthogonal hf field of 
a different frequency. The main results of this work and comparison to previous studies are 
presented in Section V. 



II. General Relations 



There are several prominent sources of nonlinear electron response in superlattices. Elec- 
tron dynamics in narrow minibands features phenomena manifested as Bloch oscillations, 
static electron localization and the Wannier-Stark ladder [|]|| in a static electric field; also, 
in a harmonic electric field there are nonlinear oscillations with amplitudes modulated by 
Bragg reflections and dynamic electron localization with the collapse of the electron's 
minibands [O . In the case of interminiband transitions, nonlinear phenomenology includes 



interminiband tunneling in a static electric field, with or without photon-assistance |23[ , also 
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interminiband tunneling in a harmonic electric field and Rabi oscillations [[24J]. Of course, 
there are also nonlinear electron response properties involved in the relaxation processes 
that redistribute energy supplied by the electric field among the various degrees of freedom, 



controlling the anisotropic heating of the electron gas [|14],[n|, upon which our attention is 
focused in this paper. 

We revisit the analysis of electron dynamics in a single miniband from a different per- 
spective than that of earlier studies |]6|j7|,p!3|1 . In this, we will examine electron dynamics in 
two-dimensional (// = 2) and three- dimensional (/i = 3) superlattices in the presence of an 
oscillatory electric field having 11 frequency components, 

E(t) = x a E a cos(uj a t -5 a ), li = 2, 3, (1) 

a=l 

where x a are the unit lattice vectors of the crystal, 8 a are the initial phases of the fields and 
the frequencies u a are different in general. We take the electron energy dispersion relation 
in the tight-binding approximation as: 

H A 
e(k) = £ a(k a ), £ a {k a ) = — (1 - cos(k a d)) , (2) 

a=l Z 

where A is the miniband width, e a and k a are the energy and wave number along the a-axis, 
respectively. 

Under the influence of the electric field of Eq.(l), an electron in a \x- dimensional super- 
lattice executes nonlinear oscillations (with a different period in each direction) having the 
velocity given by (no scattering) 

oo 

V a (k£\t ,t) = V m ^ J n (g a )sin n(u a t - 8 a ) + kf>d- g a sin(u a t ) , (3) 

n=— oo 

where V m = Ad/2h is the maximum electron velocity, k^ is the electron wave vector at 
initial time to, and g a = /u a = eE a d/htu a . One can obtain a similar expression for 
electron energy by integrating the relation V a (k a ) = h^ 1 de a (k a )/dk a . 

It is clear from Eq.(3) that the velocity harmonic amplitudes are oscillatory functions of 
the field amplitude. They sequentially vanish at the zeros of the J n (^Q,)-functions indepen- 



dently of the initial electron momentum. In particular, the time averaged values of electron 
velocity and energy (zero th harmonics) are given by: 

V a (k£\ t ) = V m J (g a )sin (k^d - g a sin(uj a t )^ , (4) 

and 

e^(k { °\ t ) = ^ 1 - J (g a )cos (k^d - g a sin(uj a t )) . (5) 

At the specific values of g a for which Jo(g a ) = 0, we have 

V*{kg\to) = 0, e^(k£\to) = j, (6) 

i.e. electron motion along the a-axis has no net translation independently of its initial 
momentum, and its average energy takes the value at the middle of the one-dimensional 
miniband. This phenomenon is known as dynamic electron localization. The electron has a 
discrete energy spectrum due to the finiteness of the motion and, therefore, dynamic electron 
localization corresponds to the collapse of its quasienergy minibands [TB| [described by the 
pre-collapse relation 

£ Q (k a ) = y [1 - J (g a )cos(k a d)] +n a hu on n a = 0, ±1, ±2, ... . (7) 

Depending on the values of g a , dynamic electron localization and miniband collapse can be 
one-dimensional, two-dimensional, or three-dimensional (in which case the localization and 
collapse are complete). 

The dynamical peculiarities of superlattice electrons are evident in nonlinear conduction. 
However, even a qualitative analysis requires the use of a correct model of the relaxation 
processes, which may be simplified to the specifications of a particular problem. In the 
present work we endeavor to take account of the /^-dimensionality of the electron scattering 
processes and the separation of elastic and inelastic scattering. The single r-approximation, 
useful for the one-dimensional model, is not adequate for our purposes, as discussed above. 
A two relaxation time model was proposed in Ref. [p5j as well as in Ref. || by one of the 



authors of the present work. However, this model is, in fact, one-dimensional and, moreover, 
it produces the illusion of a separation of elastic and inelastic electron scattering processes 
and associated scattering times, which we now understand to be incorrect. In some sense 
this model is even worse than usual single r-approximation (but, unfortunately, is still in 
use (Ref. [fZfJ)) because its identification of the two relaxation times from experimental data 
is erroneous. The balance equation method, developed in Ref. [[27|, would be useful for 
our goals, if it would be generalized by replacing the isotropic electron temperature by an 
anisotropic one (Ref. |15|), but this generalization presents a considerable challenge and has 
not been done yet. The approach proposed in Ref. [28[ takes into consideration the interplay 



between different degrees of freedom and separates phonon and impurity scatterings on a 
microscopic basis. However, this method is essentially single-particle in nature and it is 
primarily applicable for a superlattice with low carrier concentration when electron-electron 
scattering does not play a significant role. To overcome these limitations, we start from a 
three-dimensional model of a superlattice having the novel phenomenological collision term 
proposed in Ref. [15]. This collision term describes scattering in terms of an improved three- 



channel electron relaxation process. We employ the commonly understood [|T^,|r^] separation 
of the relaxation processes into elastic, inelastic and electron-electron with characteristic 
times specified for each of the three channels working in parallel. In the first channel 
(usually the fastest one) an electron is subject to redistribution of the additional energy 
and momentum supplied by the applied electric field among its degrees of freedom by means 
of elastic scattering during a characteristic time T\. The kinetic energy of each electron is 
conserved during this scattering to isoenergetic surfaces, but the direction of momentum is 
randomized (with consequent reduction of the drift velocity). In the longer-lasting second 
channel, the energy supplied by the external electric field is redistributed among all electrons 
due to inelastic electron-electron scattering, including Umklapp processes. As a result of 
the Umklapp processes the Fermi distribution becomes undrifted |T7| and, furthermore, the 



redistribution of energy establishes an effective electron temperature T e during a time r ee 
by electron-electron scattering. The total energy of all electrons is conserved during the 



redistribution process, in spite of their momentum relaxation. Finally, in the third channel, 
electrons transmit energy to the lattice over a time r e and their distribution relaxes to a 
Fermi function at the lattice temperature To. 

The Boltzmann equation with this model collision term has the form: 
df(k,t) , eE(f)3/(M) fdf\ 



dt h dk \dt 

where 



(8) 



coll 



(8J\ f(k,t)-f s (e,t) /(k,t) -f (e, r„) /(k,t)-/ (£,T ) 

UL~ * + ^ + * • () 

with the isoenergetic distribution function, fs(e,t), expressed as an integral average over 
the equipotential surface S e , 

/*(*,*)= , 'J'^' , (10) 



and 

(5) = (e) a = (e) e , (11) 

where we use the notation e = e(k) and /(k, t) is the nonequilibrium distribution function 
driven by the electric field. f (s,T e ) is the equilibrium Fermi distribution at the elevated 
electron temperature T e and fo(e, T ) is the equilibrium Fermi distribution at lattice tempera- 
ture T . (e), (s)s, (s) e , and (e) are the energies averaged over the corresponding distribution 
functions. The effective electron temperature T e is determined by Eq.(ll). It is important 
to note that the anisotropic function fg(e,t) plays the role of the "isotropic" distribution 
function of Ref. fL6| , but not the symmetric one, i.e. fs(s,t) 7^ (/(k, £) + /(— k, t))/2. A 



symmetric form of fs{e,t) would occur is there were no energy and momentum redistribu- 
tion among all degrees of freedom and it would correspond to the one-dimensional model of 



a superlattice used in Refs. H.|25]. 



Let us summarize our generalizations and simplifications of the three relaxation processes: 
(1) As in Ref. the quasielastic electron scattering leading to the " isotropization" of 



the electron distribution function is considered to be dominant. In our case, the electrons 
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are distributed onto corresponding (nonspherical!) isoenergetic surfaces and the anisotropic 
function f s (e,t) plays the role of the "isotropic" distribution function of Ref. flUfl . 

(2) The Brillouin minizones are narrow for superlattices, lending importance to Umk- 
lapp processes and the establishment of an undrifted Fermi distribution with an effective 
temperature T e W7\ . This underscores the difference of our present description from that of 



Ref. [|7|. 



(3) The dynamical development of the deviation of the electron distribution function 
from the "isotropic" one is described by 



|(/(M)-/.(M))„=- /fef) ; / - (M) . (12) 



i.e. by the effective relaxation time r, which is, in general, dependent not only on energy but 
also on the electron momentum direction. This relation is the same as the one commonly 
used for the first term of the distribution function expansion in Legendre polynomials [JTGj . 



Thus, in this respect, our description is the same as that of Ref. [16|] up to this point. In 
accordance with our classification of the three relaxation processes, we have 

I = I + ^ + I. (13) 

(4) We describe the collision dynamics of the "isotropic" distribution function f s (e,t) 
approximately by two relaxation times, r ee and r e , i.e. by the relation 

fdf s (e, t) \ f s (e, t) - / (e, T e ) /„(e, t) - /„(e, T ) 



Ot J sf T ee T s 



(14) 



While this is not an exact equation for f s (e,t), it is acceptable for a qualitative description. 
Moreover, our interest is not in the details of the distribution functions, but only in the cur- 
rent and in the average energy. Furthermore, it is easier to incorporate necessary corrections 
directly in balance equations to be derived below. 

The balance equations can be obtained from Eqs.(8-ll). For sake of simplicity we take the 
relaxation times tl, r ee , and r £ to be energy and momentum independent. Multiplying Eq.(8) 
sequentially by de a (k a )/dk a and by e a (k a ) and integrating over first Brillouin minizone, we 

9 



obtain the following equations for the current components (j a (t) = neh 1 (de a (k a )/dk c 
and the average energies (e Q ): 

- ne 2 (m-\e))E a (t) = (15) 



and 



* {£a) - l -E a {t)m = -M^i£k _ M^ifk _ ^>-^>o . (16) 

n n r ee r e 



E (Os = E (^)e = (e>, (17) 



where 



n is the three-dimensional electron density and r~ x = rf 1 + r" 1 + r^ 1 is the overall inverse 
electron relaxation time. 

The balance equations (Eqs. (15-17)) are valid both for homogeneous semiconductors 
(without superlattice) and for semiconductor superlattices of any dimensions. All terms in 
these equations have clear physical meaning and allow some generalizations, for example, 
the replacement of scalar relaxation times by a relaxation tensor (for symmetric structures 
we consider it unnecessary). It should be noted that the relaxation times, r p and r e , can 
be taken from independent calculations using the actual scattering mechanisms. Such a 
calculation for one-dimensional GaAs-based superlattices was done in Ref. |29|| , where, in 



particular, it was shown that the relaxation times, t p and r £ , can be taken to be independent 
of energy, if the miniband width is less than optical phonon energy. 

In general, the set of equations Eqs. (15-17) is not closed because of coupling to higher 
order moments of the distribution function /(k, t). However, there are simplifications for 
one-, two- and three-dimensional superlattices with a sinusoidal dispersion law, Eq.(2). In 
these cases, symmetry dictates that 

(e a ) S = (e a ) e = 7 < £ > ' ( 19 ) 
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Furthermore, the sinusoidal dispersion law provides the linear proportionality between the 
effective electron mass and its energy: 

(m-'ie)) = (y - (e a >) J- (20) 

Accordingly, Eqs.(15,16), taken jointly with Eqs. (19,20), form a closed set of equations. 
It is convenient to write this set in complex form, introducing a dimensionless complex 
yU-component "vector" with components defined by 

where jo a = (end /h) (A/2 — (e a )o). The balance equation for equivalent to 

Eqs. (15, 16, 19, 20), is given by 

+ (V 1 + «l a {t))* a (t) = r- 1 + -(r- 1 - r- 1 ) £ Re&p(t), (22) 

where £l a (t) = edE a (t)/h. The last term on the right side of Eq.(22) describes the redis- 
tribution of electron energy and momentum among all degrees of freedom. This feature is 
absent in the single relaxation time description and is of crucial importance for our present 
considerations. For a one-dimensional superlattice, (fi = 1), Eqs. (22) are identical to the 
balance equations obtained in Ref. ||. It should be noted that the only significant feature 
of the distribution functions, fs(e,t) and fo(e,T e ), is that given by Eq.(19) in regard to 
the derivation of Eq.(22) for two- and three-dimensional superlattices, because of the high 
symmetry of the electron dispersion relation. The specific forms of fs(s, t) and fo(e, T e ), be- 
yond Eq.(19), are not pertinent. A further consequence of this symmetry is that the inverse 
relaxation times rf 1 and r^ 1 are involved in Eq.(22) only in the form of their sum. For bulk 
semiconductors having a one-dimensional superlattice with free motion in the lateral plane, 
one obtains an integro-differential equation (allowing only numerical solution) instead of 
Eq.(22). Such an integro-differential equation was analyzed for the case of a static electric 



field in Ref. (1T51 in detail. 
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It is easily shown that the complex "vector" with components & a (t) introduced above is 
constituted by the first Fourier-components of distribution function f(k, t). Using periodicity 
in k-space, this distribution function may be expanded in a /i- dimensional Fourier series: 

/(k, t)=J2 F u exp{iv ■ kd}<$> v {t), (23) 

with Fourier coefficients given by 

F v = i^X J Bz d»kf (k)exp{-iis ■ kd}, (24) 



where v — > [v\) for one-dimensional, v — > (1/1,1*2) for two-dimensional, v — > (^1,^2,^3) for 
three-dimensional superlattices, and the integration is taken over the first Brillouin zone 
(BZ). Only the first momentum harmonics (in any direction) of the distribution function 

= / d^kf{k,t)exp{-ik Q d} 
a{) fd»kf (k)exp{-ih a d} 1 j 

contribute to the current density j(t) and the electron energy (e a ) for the miniband case of 
Eq.(2). One can see that they are the same as those given by Eq.(21). 

For arbitrary time-dependence of the electric field E(t), it is useful to write <3> Q (t) in a 
form that is convenient for the representation of Bloch oscillations, as 

$«(*) = a a (t)V a (t), (26) 

where 

* Q (0 =exp^-i J^VL^dt^ (27) 

is a solution of the homogeneous counterpart of Eq.(22) associated with Bloch oscillations 
in the absence of scattering terms. This function describes the dynamic modulation of the 
electron distribution function by the applied electric field without scattering. In particular, 
for a simple harmonic field, E a (t) = E^cosut, we have 

00 

* a (t)=exp{-ig a simjt}= ]T J v (g a )e- ivut , (28) 
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where g a = /uj is the projection of the dimensionless field amplitude onto the 
The functions a a (t) represent dissipative processes, describing changes in the amplitude 
of $ Q (t) (deviation from the Bloch oscillation solution) due to scattering. They obey the 
following equations, based on Eq.(22), 



a a (t) + r^a a (t) = r^Kit) + - (V 1 " ^ ( E a P*P I *«(*)■ 



(29) 



In the absence of scattering, a a (t) = 1. The transfer from the description in terms of the 
functions $ Q (t) to a description in terms of functions a a (t) corresponds to a transforma- 
tion to a new system of coordinates, K , oscillating in momentum space together with the 
unscattered electron. In the system K , each electron is at a fixed point k . Only scat- 
tering changes the distribution of these points. In the case of rare collisions {ujt » 1), 
such changes are small during the period of the applied field, but they can accumulate and 
become important over the time of a few collisions. Otherwise, the equilibrium distribution 
functions (to which the electrons relax) in the system K are modulated by the field and 
become rapid functions of time, as 

/o(e, T ) - f (h) = f (k -- d f Q fifa)^) . (30) 

This feature is embodied in the structure of Eq.(29): the "scattering-out" term (second term 
on the left side) has the usual relaxation form with the overall inverse relaxation time due 
to all scattering mechanisms, and the " scattering-in" term (first term on the right side) is 
the dynamically modulated equilibrium distribution function with inverse relaxation time 
tT 1 . The last term on the right side of Eq.(29), describing the redistribution of energy and 
momentum over the various degrees of freedom, is modulated by the electric field twofold: 
once in connection with the transformation to the system Kq, and, secondly, because the 
corresponding equilibrium distribution functions are determined by the average normalized 
electron energy (Eqs.(21),(26)) which, in turn, depends on the field and time. This average 
energy is involved in the last term through the relation 

Re (f a^ p ) = Re (f = E ^ " ~ ^ " 



p=l ' ) \I3=1 J /3=1 A / 2 ~ ( £ P)o 
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III. Self-induced transparency 



In this section, we analyze superlattice response to a high frequency {ujt p » 1) 
monochromatic field directed along the Xi-axis. Considering a a to be slowly varying func- 
tions of time, we average Eq.(29) over an interval At given by 27t/uj << At << t p , obtaining 
an algebraic system of equations, for which the stationary solution is 

„ = yjp(gi) 

1 nii+iv-m-Bigj-^-m-xyw 1 ' 



and 



n x f iB(g 1 ) + (l-\)Ji(g 1 ) 

a ~ [1 + ~ 1)A] • Sfo) - (/, - 1)(1 - A) 2 J^(g-i) ' a * A ' ^ 



where A = r p /r £ and 



£(«?) = 1 - L-^ [1 + J (2ff)] • (33) 



In the derivation of Eqs.(31),(32), we used the relations 



*«(*) = Jo(ffa), (^e(M/ Q (t)) 2 = - [1 + J (2^)] , (34) 

where the overhead bar indicates averaging over the period of the impressed electric field. 

According to Eqs. (21), (26), (28), (31) and (32), the current, energy dissipation rate, Q, 
and ratio between the transverse and longitudinal electron heating, 5, are given by 

h = hi 3m = aisin(gisin(ut)) + O (—') ; (35) 

\UJT / 



T~e T E + {fl-l)T p 

and 
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^ (/i - l)((e tt ) - (e a ) ) = (l-A)(/i-l) 
(ei)-(ei>o 1 + A(//-1) ' 



Figure 1 depicts the function ai involved in the expression for the current, and also shows 
Q/Qo, as functions of g for a three-dimensional superlattice (/i = 3) with A = 1, 0.1, 0.01. In 
the case t p = r e (A = 1), we obtain a\ = Jo(fi'i) an d expressions (35) and (36) coincide with 
those obtained for the one- dimensional model of a superlattice ||. With increasing t p while 
holding t £ fixed, the current in the superlattice increases, the dissipation decreases, and the 
amplitude of their modulation increases. As in the case of the usual single r-approximation 
|7|-||, the three-dimensional superlattice current vanishes at fields such that Jo(gi) = (with 
an accuracy of (cuTp) -1 ), where dissipation is maximal. When elastic collisions dominate over 
inelastic scattering (r p << r e ) there is a rapid redistribution of energy among the electron 
degrees of freedom and, consequently, the current decreases sharply for g\ > 1. This is 
caused by a strong expansion of the distribution function in momentum space, due to both 
longitudinal and transverse heating. 

To examine the peculiarities of self-induced transparency in the case of a two-dimensional 
superlattice, we consider an electric field of the form, 

E = (£1X1 + E 2 x 2 )cos(ujt). (38) 

In a manner similar to that of the foregoing analysis, we obtain the components of the 
dissipative function as, 

= 4A (J G?i, 2 )£(ffi,2) + (1 - A) Mgi, 2 )(Jo(9i + gg) + M9i ~ 92))) ,„ q , 
l6B(g 1 )B(g 2 )-(l-\)%Ug 1 + g 2 ) + J (g 1 -g 2 )y ' 1 ] 

To obtain Eq.(39), we used the relations (34) and the expression 

Re^JJjRe^tF) = ^ (J ( gi + g 2 ) + J ( gi - g 2 )) . (40) 



The corresponding current and energy dissipation rate are given by 

j l2 (t) = a 1>2 sin(g 1>2 sinujt), (41) 
15 



and 

Q = n{A ~ {£)o) (l - ifa J ( 9l ) + a 2 J (g 2 )^j . (42) 



It is clear from Eqs.(39),(41) that the current along each superlattice axis depends on 
all field projections, in contrast to the results of Ref. [ 30fl . In general, the current of each 



excited harmonic is not parallel or antiparallel to the resultant electric field and has its own 



elliptic polarization. There are exceptions if the electric field is in the lattice directions [10] 



or [11], in which case the field and excited current are parallel. 

Furthermore, it is evident from Eqs.(39),(41), that both dynamic localization and self- 
induced transparency can be either one- dimensional or complete for two-dimensional su- 
perlattices. In Figure 2 we show the locii of one- dimensional self-induced transparencies 
occuring in directions xi (horizontal curves) and X2 (vertical curves), respectively, as func- 
tions of gi and g 2 ■ At fields corresponding to these curves, the current components vanish 
CM*) ~ and j 2 (t) 0, respectively). These dependencies exhibit oscillations around the 
line Jo (01,2) = (see Figure 2 Insert), which is the condition for dynamic localization to 
occur. The average energies and the dissipation rate for these values of the fields are given 
by 

A 7 r_ A f A i v\ 4AJ 2 (£ 2{1) 



<eim> = _ M = T -^-W.J 4+(1 _ A) — (W . (43) 



and 



Q = "< A -< £ >°> 1 2A f^>> ,1, (44) 

' 4 + (l-A)(l + J„(2 92(1) :^ 



where the first of the subscripted indices is related to the horizontal curves and the second 
subscripted index (in parentheses) is related to the vertical curves. At points of intersection, 
determined by the condition Jo(gi) = Jo(92) = 0, there is complete self-induced trans- 
parency. In such cases, the average electron energies and the dissipation rate are maximal 
and are given by 
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(ei> = (S2) 



A 

2"' 



(45) 



and 



Q = 



n(A- (e) ) 



(46) 



The dissipation rate in Eq.(46) is larger than that in Eq.(44), and is also larger than the 
maximum dissipation rate for the fields oriented strictly along the crystal axes, the latter 
being given by Eq.(36) as 



It is apparent that, within the three-relaxation-time description presently under consider- 
ation, complete self-induced transparency and dynamic localization occur at the same fields, 
whereas one-dimensional self-induced transparency and dynamic localization arise at differ- 
ent fields. In two-dimensional superlattices complete dynamic localization and self-induced 
transparency occur at the discrete amplitude values and applied electric field directions 
determined by the relations 



where (p m ,n is the angle of field orientation with respect to the superlattice crystal axis and 
£ m is the m'^-order root of the zero'^-order Bessel function. This can be generalized easily 
to the three-dimensional case. 

To develop a physical understanding of one-dimensional and complete self-induced trans- 
parencies, we examine the evolution of the electron distribution in the ii'o-system. One can 
see from Eq.(29) and its following discussion that, under the influence of the field and 
scattering, the number of electrons entering the current components, Im^i^ (including re- 
distribution among them), averaged over the field period, P i;2 , is given by 



P 1>2 = T^i, 2 (t) + - (r" 1 - t- 1 ) (a h2 (ReV h2 (t)) 2 + a 2yl Re^ 1 (t)Re^ 2 (t)) . (49) 




T £ + T p 



(47) 




(48) 
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If the dynamic modulation of the equilibrium distribution function is such that ^i(t) = 
^2(t) = 0, then the average number of electrons entering the current components van- 
ishes (there is only electron redistribution among the components), and, therefore, the 
current components are eventually completely eliminated from the nonequlibrium distri- 
bution function by "scattering-out" over a time of order t p . It should be noted that for 
^i(t) = ^2(t) = 0, the set of equations (29) averaged over the field period becomes homo- 
geneous and its steady-state solution is zero. In this case electron heating is maximal due 
to complete dynamic localization. Thus, after a time of order t p , the superlattice becomes 
transparent, i.e. behaves like a dielectric having the permittivity of the crystal lattice and 
relatively small, but resonant, absorption. This is to say that we have complete self-induced 
transparency. It should be emphasized that, at arbitrary fields, the absorption rate stabilizes 
after time r £ , i.e. later than the vanishing of the current. 

If dynamic localization takes place in only one of the two crystal directions, for example, 
in [lOj, i.e. ReW^Jt) = J (gi) = 0, but Re^ 2 {t) = J (g 2 ) ^ 0, then the electrons only 
enter the current component §2(t) due to dynamic modulation of the equilibrium distribu- 
tion function. However, because of the redistribution of energy and momentum among the 
degrees of freedom, electrons also flow into <J> 1 (t). Therefore, even at Jo(<?i) = both P 1 
and ji are nonzero. The current ji vanishes only if the components of the flow P 1 caused 
(a) by direct dynamical modulation of the distribution function, and (b) by redistribution 
via scattering, compensate each other in this direction. However, this occurs at Jo(fl'i) 7^ 0, 

1. e. when dynamic localization in this direction is absent, as reflected in Eq.(32) and Figure 

2. Furthermore, it should be noted that the time-averaged one-dimensional energy is A/2 
for each electron when dynamic localization occurs, but this is valid only for both time and 
ensemble averaged energy in the case of one-dimensional self-induced transparency. 

One can easily see that if the relaxation times r v and r e are energy dependent {not con- 
stant), then the first harmonics of the distribution (23) become coupled to higher harmonics. 
In this case self-induced transparency does not occur even when there is complete dynamic 
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localization. This was demonstrated in Ref. |7[ for a one- dimensional superlattice. How- 
ever, as mentioned above, for an appropriate set of superlattice parameters, these relaxation 
times can be considered energy-independent |29j and, therefore, even the quantitative results 
discussed above have a wide range of validity. 

It should be emphasized that, although both self- induced transparency and dynamic 
localization occur in a superlattice due to the narrowness of its Brillouin minizones, the 
physical origins of these effects are completely different. Dynamic localization arises when 
the zero th harmonic of nonlinear electron oscillations, modulated by Bragg reflections, van- 
ishes. In contrast to dynamic localization, self-induced transparency is a result of the joint 
action of Bragg reflections of miniband electrons and collisions creating strongly modulated 
electron distributions in which the first harmonics are absent for discrete values of the elec- 
tric field amplitudes (in the limit r v — > oo). As a result, dynamic localization appears 
immediately after turn-on of electric field, and self-induced transparency occurs only after 
a time of order of t p . The conditions for self-induced transparency depend on the scattering 
mechanisms and, in general, it takes place even without dynamic localization. It can be 
shown that states of self-induced transparency are not stable with respect to the generation 
of static and hf fields (having frequencies not equal to u>i) ]3T|] and, therefore, experimen- 
tal studies should be performed at low electron concentrations and with the use of pulsed 
electric fields. 



IV. Current modulation by orthogonal fields 



In this section we examine superlattice behavior in the presence of a high frequency 
biharmonic electric field given by 

E = Eix. 1 cos{u)it — 5i) + E 2 ^2Cos(u2t — (50) 

In this case the field components have different frequencies and are directed along different 



crystal axes taken to be orthogonal (the case of parallel fields was analyzed in Ref. |32fl in a 
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single r-approximation). We are interested to explore the occurrence of current amplitude 
modulation by a high frequency electric field orthogonal to the current direction. This is 
determined by the redistribution of energy and momentum among the various degrees of 
freedom. We assume that the frequencies of the most important electric field harmonics are 
well separated, i.e. 

\niUi - u 2 u; 2 \t p » 1, n lj2 — 1, 2, ... . (51) 
Using Eq. (34) and the expression 



Re^ a {t)Re^ p {t) = Mg a )J (g fi ), a^(3, (52) 
we obtain the following relations for current components and averaged energies: 

ji,2 = a 1 fl(g 1 ,g 2 )sin(gi i2 sin(ui i 2t - <$i, 2 )), (53) 

and 

(£1,2) - (£1,2)0 = (y ~ (£1,2)0) (1 - a h2 (g 1: g 2 )Jo(gi, 2 )) , (54) 

where, in the case of a two-dimensional superlattice, 

„ (n , \ _ 4AJ (g 1|2 ) (gg^i) + (1 - A) Jo(ff 2 ,i)) ( ,,, 

whereas, for a three-dimensional superlattice we have 
01,2(01,02) = 

9AJ (^, 2 ) (31^,0 + (1-A)J 2 Q ?2 , 1 )) 

9(2 + \)B( 9l )B(g 2 ) - (1 - A) 2 (3J$( gi )B(g 2 ) + 3J*(g 2 )B( gi ) + (4 - A) J*( 9l ) J 2 (^ 2 )) 



(56) 



In Figure 3a we exhibit the functions a x 1 (0.05, g 2 ) = a\{g\ = 0.05, g 2 = 0)/ai(gi = 
0.05, g 2 ) and a^i, 2.405) = a 1 (g 1 ,g 2 = 0)/a 1 {g 1 ,g 2 = 2.405) for A = 0.1. The former 
function describes the modulation due to the electric field E 2 (orthogonal to E\) of current 
driven by a weak electric field E\\ whereas the latter function describes the same current 
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modulation in the presence of an arbitrary field E\ under conditions of dynamic localiza- 
tion in the direction x 2 (Jo(fl , 2) = 0). The functions a 1 (g 1 ,g 2 = 0) and a 1 (g 1 ,g 2 = 2.405), 
corresponding to the curve a^ 1 ((?i, 2.405) of Figure 3a, are shown in Figure 3b. 

One can see from Eqs.(55),(56) and Figures 3a,b, that it is impossible for the orthogonal 
field to cause a complete vanishing of polarization, and, therefore, induced transparency 
does not occur. However, the modulation of polarization by the orthogonal field can be 
significant, especially for small A. The maximum decrease of polarization comes about with 
the occurrence of one-dimensional dynamic electron localization in the transverse direction 
{Jq(92) = 0). In particular, for E\ — > 0, it decreases by the factor (1 + (/i — l)A)//iA. The 
reason for this is that, under condition of dynamic localization, electron heating is maximal 
and the distribution function widens in all directions in momentum space due to energy and 
momentum redistribution among all degrees of freedom, which always leads to decreased 
current. 

As in the case of self-induced transparency, for specific ratios of field amplitudes and 
frequencies determined by the conditions Jo(gi) = Jo(g 2 ) = 0, complete induced superlattice 
transparency takes place. In this case a two- dimensional superlattice is transparent to an 
arbitrarily polarized third weak signal with frequency u 3 (well separated from the frequencies 
riiUi ± n 2 u 2 , where ri\ and n 2 are integers). 

Similar to self-induced transparency, current modulation by an orthogonal field and 
induced transparency occur (vanish) in a time of order of r p and they become stationary in 
a time of order of r e after turn-on (turn-off) of the electric field, i.e they are displaced in 
time from dynamic localization. 

To establish the frequency limitations of the phenomena we have explored, we carried 
out a numerical analysis of Eq.(22) with finite values of r p and t £ . The results are shown in 
Fig. 4, where the dashed lines represent the above- described analytical calculations for the 
functions 1 (0.05, g 2 ) (upper curve) and a^ 1 ((?i, 2.405) (lower curve), respectively. The solid 
lines represent our numerical determinations of the amplitudes of the first current harmonics 
(thick curves) and the maximal currents (thin curves) at ur p = 1 (maintaining ujt £ » 1). 
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One can see that the dependencies presented in Figure 3(a,b) change only slightly with 
decreasing t p . 

V. Summary 

In summary, we have applied the Boltzmann equation with an improved three-relaxation- 
time collision term to the analysis of self-induced and induced transparencies in semiconduc- 
tor superlattices. The three relaxation times include (a) a time for redistribution of energy 
and momentum supplied by an electric field to a given electron among its various degrees 
of freedom, (b) a time for redistribution of energy and momentum among all electrons by 
inelastic electron-electron interactions, and (c) a time for transfer of the excess energy to the 
crystal lattice. We have performed analytical calculations for systems having high symme- 
try (for one-, two- and three-dimensional superlattices). However, the results obtained here 
are valid qualitatively for bulk semiconductors with a one-dimensional superlattice, which 
are currently available for experimentation. Furthermore, we have shown that self-induced 
transparency and dynamic localization are different phenomena with different physical ori- 
gins, displaced in time from each other, and, in general, they arise at different electric fields. 
Moreover, we have found that the redistribution of energy and momentum among the various 
degrees of freedom is of crucial importance in two-dimensional and three-dimensional super- 
lattice transport and optical properties. Transverse electron heating drastically changes the 
conditions for self-induced transparency, and this effect facilitates current modulation by an 
applied perpendicular high-frequency electric field. 
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Captions to the Figures: 



Figure 1. Normalized (a) current and (b) energy dissipation rate as functions of g = 
eEd/hu. 

Figure 2. Locii for one-dimensional self-induced transparencies in directions x x (hori- 
zontal curves) and x 2 (vertical curves). The Inset exhibits a magnification of the locus for 
one-dimensional self-induced transparency in direction x x and the locus for one-dimensional 
dynamic localization in this direction (horizontal line). 

Figure 3. (a) Current modulation by orthogonal fields as described by of 1 (0.05, g 2 ) 
and aj" 1 ((7i, 2.405) [defined in text]; (b) The functions ai(gi,g 2 = 0) and ai(gi,g 2 = 2.405) 
employed in the determination of a^ 1 ((7i, 2.405) of Fig.3a. 

Figure 4. Amplitude of j/jo as a function of field (#1(2)) • 
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